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Abstract. We argue that the textbook method for solving eigenvalue equations 
is simpler, more elegant and efficient than the Asymptotic fteration Method (AIM) 
applied in J. Phys. A 44 155205. We show that the Kratzer potential is not a realistic 
model for the vibration-rotation spectrua of diatomic molecules because it predicts 
the position of the absorption infrared bands too far from the experimental ones (at 
least for the HCl and H 2 molecules chosen as illustrative examples in that paper). 



In order to study the vibration-rotation motion of diatomic molecules in N 
dimensions Durmus [1] chose the Kratzer potential and solved the Schrodinger equation 
in hyperspherical coordinates by means of the Asymptotic Iteration Method (AIM). 
He obtained the well known results and as an illustrative and practical application of 
the model he restricted himself to the only apparently relevant case N = 3 calculating 
some vibration-rotation energies for Hi and HCl. In what follows we contrast the 
AIM derivation of the main equations with the well known and widely used textbook 
approach, and also compare the theoretical results for those molecules with experimental 
ones. 

The starting point is the Schrodinger equation 



where \i is the reduced mass of the molecule and E is the vibration-rotation energy. As 
a "realistic" model for the interaction between the nuclei the author chose the Kratzer 
potential 



where r e is the equilibrium internuclear separation and D e is the molecular dissociation 
energy (misleadingly called intermolecular separation and dissociation energy between 
diatomic molecules, respectively, by the author p]). He made a curious distinction 
between the modified Kratzer potential 77 = and the Kratzer potential rj = —D e which 
are just two alternative expressions of the same interaction with the energy origin shifted 




(1) 




(2) 
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by 77. Without a plausible justification he further argued that it is useful from a physical 

point of view to consider the general iV-dimensional case where 
N N pa 

'' = E4V = e4 (3) 

j=l j=i ux j 

Since the potential is spherically symmetric one can separate the Schrodinger 

equation ([1]) into its radial and angular parts, and by means of well-known 
transformations the former is reduced to pQ 

yF"(y) + (2 7 + N - 1 - 2f3y) F'(y) + [2k - /3(2 7 + N - 1)] F(y) = (4) 

where 

r 

y = — 
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q2 2^ e (D e -E + V ) 
P ~ h 2 



N-2 / / iV-2\ 2 
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and / = 0, 1, ... is the angular-momentum quantum number (called In-i in Ref. PQ). In 
order to solve this equation Durmus [1] applied the AIM which is an iterative approach 
that gives results for n = 0, 1, . . ., where n labels the number of iterations. By inspection 
of the particular outputs j3oi, (3u, . . . one hopefully derives the value of (5 n i for an arbitrary 
number n of iteration steps and then the allowed energy E = E n \. With some more 
ingenuity one realizes that F(y) is proportional to the confluent hypergeometric function 
n, 27 + N — 1; 2/3?/). This procedure offers little difficulty if one already knows 
the exact result beforehand which is actually the case here. Once we have the solution 
in terms of the confluent hypergeometric function we easily rewrite it in terms of the 
associated Laguerre polynomials [UE]. 

Solving equation (j3J) is a textbook problem [2] and the widely known approach is 
faster, more elegant and efficient than the AIM. If we define y = ^ then w(z) = F(j^) 
satisfies 



zw"{z) + (2 7 + N — I — z) w'(z) + 



k ( N 



w(z) = (6) 
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that is a particular case of Kummer's equation [3] 

zw"(z) + (b — z) w'(z) — aw(z) = (7) 

The Kummer's function can be expanded in a power-series 

j=0 J-°j 

where = £(£ + 1)(£ + 2)...(£ + j — 1) and £o = 1- It becomes a polynomial of degree 
n when a = — n and 6 7^ — m (m and n positive integers). On comparing equations (JH]) 
and © we directly obtain 

and [3J 

F n {y) = M(-n, 2 7 + N - 1, = iFi(-n, 2 7 + iV - 1; 2/3y) (10) 

It is clear that one can derive the solution to equation (j4j) directly from comparison of 
the appropriate equations in a way that makes the AIM utterly unnecessary. 

Durmus pQ did not show any physical application of the iV-dimensional model 
except for the obvious case N = 3. Most curiously he showed results in the form of a 
table and figure for both r/ = and rj = —D e . Apparently, he did not realize that both 
spectra are the same but for an energy shift that is irrelevant from a physical point of 
view. Suffice to say that the physical observables are not affected by this shift. 

Another curious fact is that in Table 2 the author only showed the energies for 

n > I = 0, 1, In the case of a diatomic molecule, to which the model is supposed to 

apply, n is the vibrational quantum number v and / is the rotational quantum number 
J so that such selection of quantum numbers is of scarce utility from a spectroscopic 
point of view (one would expect increasing J for a given v). 

Durmus [1] argued that the Kratzer potential provides a realistic description 
of molecular vibrations, but it is far from true as discussed by Plfva [4J who 
stated that "However, in its basic form, this function provides only a rather crude 
approximation for the molecular potential, and for this reason it has not been popular 
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with spectroscopists" . In the first place, the Kratzer potential supports an infinite 
number of vibration-rotation levels which is not the case of actual diatomic molecules. 
In addition to it, it does not describe the spectrum correctly as we will show in what 
follows. From the vibration-rotation energies written in the usual spectroscopic way 

E v ,j = (ll) 



we obtain the spectral lines in wavenumber units 

~v = E "> J '- Ev > J (12) 
he 

according to the selection rules Av = ±1 and A J = ±1. [5] The results in Durmus' 

Table 2 do not even allow us to obtain the P and R branches of the fundamental 

absorption band (v — 0, v' — 1, J' — J ± 1). [5J The first two entries for HCl in 

that table predict the center of this band to appear at E\ q — E = 0.1482 eV or 

1195 cm -1 . However, it is well known that the center of the fundamental band is located 

at v — 2886 cm -1 , [5 J more than twice the value given by Durmus' energies. It is not 

better for H2 (theoretical= 2715.7 cm -1 , experimental 4160.2 cm -1 ) as expected from 

a model that no spectroscopist would take seriously. 

Summarizing the main conclusions of this comment we may say that the application 

of the AIM for solving the Schrodinger equation with the Kratzer potential is far from 

being the best strategy. One obtains the same expressions in a more direct, easy and 

elegant way by means of the standard textbook method. |2j In addition to it, since 

Durmus did not show any plausible physical application of the model to space dimensions 

other than N = 3 and the results for ordinary diatomic molecules are extremely poor 

we conclude that the model is of scarce physical utility. 
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